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Abstract 

We address the issue whether ARPES measurements of the spectral function 

Ak{ijo) near the Fermi surface in the normal state of near optimally doped 

cuprates can distinguish between the marginal Fermi liquid scenario and the 

spin-fluctuation scenario. We argue that the data for momenta near the Fermi 

surface are equally well described by both theories, but this agreement is 

nearly meaningless as in both cases one has to add to a large constant 

of yet unknown origin. We show that the data can be well fitted by keeping 

only this constant term in the self-energy. To distinguish between the two 

scenarios, one has to analyze the data away from the Fermi surface, when the 

intrinsic piece in 5] (a;) becomes dominant. 
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The observed discrepancies between the normal state properties of the cuprates and the 
predictions of Fermi hquid theory continue to attract a lot of attention from the condensed- 
matter community. In a Fermi liquid, the quasiparticle damping near the Fermi surface 
behaves as lo^ or {irTy, whichever is larger. This gives rise to the tu^ (T^) behavior of the 
width of the peak in the quasiparticle spectral function Ak{Q) = {l/7i)ImG{k, fl) measured 
in the angle-resolved photoemission (ARPES) experiments. 

The data on optimally doped cuprates show several features which are in disagreement 
with Fermi liquid theory: T,"{uj,T) is large and does not display a Fermi liquid form for all 
T > Tc and for all frequencies. Instead, over a wide range of u and T, T,"{uj, T) resembles a 
linear function of both arguments [0-^. As a result T,"{uj, T) does not become progressively 
smaller than u with decreasing uj, and the spectral function does not undergo a drastic 
sharpening at the smallest energies, as it does in a Fermi liquid. Whether this effect is 
quantitative or qualitative is still a matter of debate as in the normal state experiments, 
T remains finite and larger than T^. The data for overdoped materials do show a stronger 
sharpening of the spectral function near kp indicating that these materials are likely Fermi 
liquids From this perspective, the unusual behavior of optimally doped cuprates indicates 
that the upper energy cutoff for the Fermi liquid behavior rapidly reduces down from the 
Fermi energy. This is a typical situation for system near a quantum phase transition. 

In this communication, we compare two scenarios for the non-Fermi liquid behavior, 
both are based on the idea of quantum criticality. The first is the marginal Fermi liquid 
(MFL) scenario developed in Ref 0. This scenario assumes that there exists a quantum- 
critical point (QCP) separating the Fermi hquid and the pseudogap regimes (Fig. la). 
Superconductivity develops in the shaded region about the QCP. In the quantum-critical 
regime, the system is assumed to display a marginal Fermi liquid behavior with 

E'^^^u) = -A„^4a;log(-) - z^Vuj^ + 7r2T2] (1) 

where Xmfi and Uq are the two adjustable parameters. To satisfy Kramers-Kronig relations, 
uo should possess a nonsingular dependence on (vrT/cu)^, but this dependence is usually 
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neglected. The self energy is defined such that G(k, cu) = — ek + S(k, c<j)). We plot 
S''"'^'(u;) in Fig. lb and c. This marginal behavior naturally emerges near the critical points 
in 3D, however obtaining this behavior in 2D systems is still problematic. To do so, it 
is possible that special requirements such as, e.g., singular momentum dependence of the 
low-energy bosonic mode which mediates the fermionic interaction [0 must be satisfied. 

Another scenario is based on the idea that the dominant interaction in the cuprates is 
between fermions and their low-energy collective spin excitations. This physics is described 
by the spin-fermion (SF) model. In this scenario, the non-Fermi liquid behavior in the 

normal state is also associated with the closeness to a critical point, but this point now sepa- 
rates paramagnetic and antiferromagnetically ordered phases. The pseudogap phase emerges 
in the SF model as a dome around this critical point, and optimal doping roughly corre- 
sponds to a situation when this dome intersects the crossover line separating Fermi liquid 
and quantum-critical regions (see Fig. Id). Superconductivity results in the shaded region 
where the pseudogap and Fermi liquid regimes intersect, and the normal state non- Fermi 
liquid behavior occurs in the quantum critical regime that near optimal doping stretches 
down to almost T^. The quantum-critical behavior in the SF model has been analyzed in 
detail. 0] The model has two parameters; ujgf which is the typical frequency for the relax- 
ational spin dynamics and also turns out to be the upper cutoff of the Fermi liquid behavior, 
and the dimensionless coupling A^/. Both A^/ and uj~j increase as the system approaches 
a magnetic instability, but the product uo = AX'j.j-uJsf remains finite and serves as the upper 
cutoff for the quantum-critical behavior. The separation into Fermi liquid and quantum 
critical regions makes sense when A^/- > 1 or equivalently when uj^f <^ u. In general, Xsf 
and Usf (but not Co) depend on the momenta along the Fermi surface, but this dependence 
is rather mild as we will show. 

The microscopic nature of the SF model implies that the self-energy is totally determined 
and has to be obtained in explicit calculations. At T = 0, the Eliashb erg-type calculations 
yield 
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K' i'^) = Kf (2) 

At small frequencies, uj <^ Usf, {uj) = Xgf {uj + iu!\u!\ /AuJsf) and the system dis- 
plays a Fermi liquid behavior. At frequencies above 8 — lOciJsf, the fermionic self-energy 
displays a quantum-critical behavior, in which u is the only energy scale: J^'lf {u) = 
gi7r/4gjgj-^^ j^^|^|-ji/2^ intermediate frequencies, T^lf (lu) interpolates between the two 
limits. Although there is no intermediate asymptotic (i..e., the crossover happens at en- 
ergies 0{ujsf)), the intermediate regime is rather wide, and in between uJsf/2 and at least 
6 — Susf, (uj) resembles a linear function of u. The finite temperature result is more 
complicated but posesses the same basic characteristics as the zero temperature result. 
We plot E'J (to) at T = lOOK in Fig. Id. 

In light of the strong similarities between the MFL and SF fermionic self energies, it 
is important to ask which experiments allow us to distinguish between the two scenarios. 
We focus in this communication on energy distribution curve (EDC) ARPES experiments 
performed for momenta near the Fermi surface. Abrahams and Varma |T^ have argued that 



such ARPES experiments help put strong constraints on the form of the pairing interaction 
and can rule out many of the scenarios for the non- Fermi liquid behavior in the normal state. 
In particular, they argued (i) that the MFL form of ^{uj) nicely fits the experimental data, 
and (ii) the data along the Fermi surface can be fitted assuming that the MFL coupling 
Xmfi is almost independent of the position on the Fermi surface. They argued that this 
rules out a magnetic scenario since the bosonic interaction is peaked aX Q = (tt, tt) and 
hence the coupling constant should be momentum dependent with the maximum at hot 
spots (points at the Fermi surface separated by Q) in contradiction with (ii). We, on the 
contrary, demonstrated previously [|11|] that the photoemission data near (0, tt) and along 
the zone diagonal can be well fitted by the spin- fluctuation self-energy, Eq. (Q), with A^/ 
varying along the Fermi surface. 



We re-examine this issue using recent photoemission results from the Argonne group. 



We will argue that EDC data for k kp can be reasonably well fitted by both the MFL and 



SF theories, but the way in which the fitting works poses a serious question for both scenarios. 
Namely, in both fits, one has to add to the self-energy a large imaginary constant i'-){kp) of 
yet unknown origin. We show that the fits with = i'-flkp) with no frequency dependent 

terms added turn out to be of the same quality as the fits with Zl{"''(ti;) = T,k{uj) + i'y where 
T,k{uj) has either MFL or SF form. From this perspective, the ARPES measurements near 
the Fermi surface do not allow the ruling out of either scenario. The situation may be 
different at large deviations from kp when the maximum in the photoemission intensity 
moves to high frequencies, and Sfc(co') (which increases with u) becomes larger than 7. In 
this case, one could possibly distinguish between MFL and SF scenarios by examining both 
the frequency dependence of and the sign and magnitude of S'(to'). 

Before we proceed with the fitting, we briefly review the experimental 
setup [P,pt^ p!2| , p!4| , p!3| , p!5| . Two types of ARPES experiments are currently available. The 
EDC (energy distribution curve) experiments measure fermionic Ik{uj) = Ak{uj)nF{uj) as a 
function of frequency at a given k. The MDC (momentum distribution curve) experiments 
measure Iki^^) as a function of k transverse to the Fermi surface at a given frequency. In 
general, near the Fermi surface, 

. = 1 

' (^ + SU^) - Vk{k - k,)f + (S'^(^))2 ^ > 

There are both experimental and theoretical indications that in the cuprates, the 
fermionic self-energy is almost independent of the momentum component transverse to the 
Fermi surface although it strongly depends on frequency, and also possesses some depen- 
dence on the momentum along the Fermi surface. In this situation, the self-energy in the 
MDC measurements remains almost a constant, and the photoemission intensity evolves 
with k as 

Ik^'^'i^) (^^(A;_A;*(cu)))2 + (S")2 ^'^^ 
where k*{uj) = kp — {uj + Tj'{uj))/vf- We see that Ik has a Lorentzian form with the 
maximum aX k = k* and FWHM \k — k*\ = Tj"{uj)/vp. In other words, the MDC measure- 
ments allow one to directly infer the frequency dependence of the fermionic self-energy. If 



vp were known, the MDC measurements would also yield the magnitude of However 
to obtain f^r is itself a problem as the MDC dispersion extracted from the peak position 
at low enough oj is uj{k) = v*p{k — kp) where v*p = vf{1 + dT,' (uj) / duj) , hence measuring 
the dispersion at low frequencies one obtains a renormalized Vp, not a bare vf One 
could hope to extract vp from the dispersion at higher frequencies, where u > S'(c(j), and 
Vp approaches vp- However, it is not quite clear whether T,'{uj) can be fully neglected at 
energies where one can restrict with the linearized dispersion near the Fermi surface. 
In the EDC experiments, one measures at small enough frequencies 

^'''''^"'^ {u;-u{kwi%"{uj)/v*py 

In a Fermi liquid S"(co') is small and near u = uj{k) can be approximated by a constant - its 
value at uj = uj{k). I^^'^{uj) then has a Lorentzian form with FWHM Tj"{uj)/v*p. As v*p is 
directly measured from the dispersion, one could straightforwardly extract from the 

data. In the cuprates, however, the situation is more complex as (i) typical frequencies are 
not small, and hence v*p itself depends on frequency, and (ii) S"(c<j) is not small, and its 
frequency dependence matters. Both effects give rise to a non-Lorentzian form of I^^^{uj), 
and the extraction of S"(co') becomes problematic. 

In Figs @, § and |^ we present the data for /fe(c(j) obtained from EDC experiments by 
Kaminski et al |jl6| taken at lOOfT from near optimally doped Bi2212 with a of 85i^'. 
We have chosen lineshapes for two points near the Fermi surface, one is close to the zone 
diagonal {ky = 0.547r/a) , and another near a hot spot {ky = O.Slvr/a). We see that the form 
of I^^'"{uj) is clearly different from a Lorentzian, even when the Fermi function is taken into 
account. In this situation, it is more reasonable to fit the whole I^^'-''{uj) rather than the 
FWHM. This is what we will do. 

In Fig.0 we present the theoretical fits using the MFL form of the self-energy. The fit in 
Fig. ^ is obtained using uq = 500meV. For this c^o, S'(ci;) is small at frequencies of few 
hundred meV (see Fig. |I|) and to first approximation can be neglected. The fit in Fig. |^ is 
obtained using a deliberately large ujq = 5eV, in which case the real part of the self-energy 
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is strong. We see that both fits are rather good. The fitting parameters in Fig ^ a (b) are 
^mfi = 0.17(0.17), 7 = 100(133)me\^ and = —26{—30)meV for k along the zone diagonal, 
and Xmfi = 0.15(0.15), 7 = 171(225)mel^, and = — 43(— 47)mey for k close to a hot spot. 
To account for the flattening of the measured I^^^{uj) at the highest frequencies, we also 
had to add a background contribution to I^^'-^{uj) which we have chosen, following Ref. |Tp 



to be just proportional to the Fermi function: I^^'^ioj) — > I^^^iuj) + {Snpioj). We found 
13 = 1.57(1.64) along zone diagonal and /? = 2.22(2.26) near the hot spots. 

We next consider the fits by the SF model. As stated previously, the two input pa- 
rameters in the model, Ugj and Xgf generally depend on the momentum component along 
the Fermi surface. Near the hot spots, this dependence is universal and is given by 
Xsf{k) = Xsf/Zl^'^, oOsfik) = uJsfZk where = 1 + {6k^y, and ^ and 6k are the spin 
correlation length and the momentum deviation from a hot spot along the Fermi surface, re- 
spectively 0. The largest 6k^ is for k along the zone diagonal. At optimal doping, ARPES 
measurements of the Fermi surface yield 6k^^^ ~ 0.27r/a ^ 0.6/a [|1],|2[. The correlation 
length extracted from the NMR measurements is ^ ~ 1 — 2a. Then Xsf is reduced by 



at most a factor of 2 as one moves from hot spots to zone diagonal. For Usf, the use of the 
expansion formula near hot spots would yield a bigger increase, but its actual variation is 
much smaller as the theoretical Ust oc sin (po, where 0o is the angle between Fermi velocities 
at k and k + Q, which tends to tt as k approaches the zone diagonal. p| In Fig. ^we present 
the results of a computation of ujsf{k) along the Fermi surface which includes both of the 
above effects. The increase is only a factor of 2.4, consistent with the fact that ^^sf{kdiag) 
extracted from fitting v*p by the SF self-energy yields ojsfik^iag) ~ 20 — 2'bmeV while 
the fits to NMR yield ujsf ~ IbmeV for fermions near hot spots |1^. For simplicity, in the 



SF fits, we will consider to be independent of k and set = 15meV in accordance 
with NMR. We, however, have verified that variations in have little effect on ARPES 
lineshapes. 

The fits using the SF form of the self-energy are presented in Fig ^ The two fits are 
obtained using momentum independent and momentum-dependent coupling Xsf{k), respec- 



tively. We see that the fits are of the same quahty as those in Fig. |^. The fitting parameters 
in Fig|a (b) are A,/ = 0.84(1.0), 7 = 110{lU)meV, p = 1.68(1.63) and = -32(-30)mel^ 
for k along the zone diagonal, and A,/ = 0.84(1.5), 7 = 180(193)meF, /3 = 2.31(2.28) and 
efc = —52{—A5)meV for k close to a hot spot. We verified that making A^/ a bit larger or 
allowing it to vary more over the Fermi surface does not change the accuracy of the fits. 

Since the frequency dependence of the self energy appears to make little difference in 
the quality of the fits at the Fermi surface, we now attempt to fit the data by S(c(j) = i'y{k) 
without extra MFL or SF frequency dependence. The results of these fits are presented in 
Fig.^. The fitting parameters are 7 = 94meV, (3 = 1.8 and = —AAmeV for k along the 
zone diagonal, and 7 = 160meV^, f3 = 2.4 and = —12meV for k close to a hot spot, is 
much smaller than in the other two fits, due to the lack of a real part of the self energy in 
these fits. 

All three fitting procedures work surprisingly well. The facts that the data are not 
obtained exactly at the Fermi surface (i.e., ek 7^ 0) and that one has to add the background do 
not differentiate between the fits and hence are not that relevant, particularly as ~ 50meV 
corresponds to a very small \k — kpl ~ ek/vp ~ 0.02(7r/a). The two physically relevant 
parameters are the coupling constant A and the extra '^{k). By adjusting 7(A;) one can fit 
the data quite well by both MFL and SF forms, and also by just S" = i'y{k). This indicates 
that one cannot differentiate between theoretical scenarios by analyzing EDC data obtained 
near the Fermi surface. It is possible that the intrinsic piece of the self energy may be 
extracted by fitting either EDC or MDC at deviations from the Fermi surface of at least 
a few hundred meV. Finally, the large '^{k) in the fits poses a problem for both MFL and 
SF scenarios. The origin of the 7(A;) is currently unknown and its understanding is clearly 
called for. 
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FIG. 1. Comparison of the marginal fermi liquid and spin fermion scenarios. Figures a and d 
show the proposed phase diagrams for the MFL (a) and SF (d) pictures. Descriptions of the phase 
diagrams are in the text. Figures b,c and e,f compare the MFL and SF self energies at T = lOOK. 
Figures b and c show the real and imaginary parts of the MFL self energy respectively. \mfi =0.1 
in both figures. In figure b, S' is given for both a;o = 500meF and = SeF. Figures e and f show 
the real and imaginary self energies for the spin fermion model. A^/ = 1.0 and a;^^ = 15meV. S" 
is seen to be roughly linear over a wide frequency range. 
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FIG. 2. Fits of the measured EDO dispersion by the marginal Fermi liquid form of the 

self-energy. The intensity /(a;) is in arbitary units. The two fits have almost the same cou- 
pling Xmfi but different loq and different 7. In Fig. a, loq = 500meV in which case the real part of 
the self-energy is small for the relevant frequencies. In Fig. b, ujq = 5eV, and S' is relevant. Both 
fits work very well. The parameters are presented in the text. 
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FIG. 3. The behavior of the spin-fluctuation frequency LOgf along the Fermi surface. This 
frequency is the upper cutoff for the Fermi liquid behavior near an antiferromagnetic instability. 
We used e{kx, ky) = — 2t(cos kx + cos ky) + At' cos kx cos ky — n with t = 250meV, t' = — 0.36t, and 
/J, = —l.lt chosen as being representative of optimal doping. The variation of tOgf at a deviation 
from hot spots is determined by a competition between the increase due to a reduction of the 
magnetically mediated interaction, and the decrease due to the "nesting tendency" - the increase 
towards tt of the angle between the Fermi velocities at k^? and k^T^ -|- Q. 
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FIG. 4. Fits of the measured EDO dispersion by the spin-fermion model with ojgf = 15meV. 

(a) a fit using momentum independent coupling Xsfik) = 0.84 for both the diagonal and the hot 

spot, (b) a fit using momentum dependent coupling Xgf = 1.0 along the diagonal and Xgf = 1.5 

near the hot spots. Other parameters are presented in the text. 
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FIG. 5. Fits of the measured EDC dispersion by the self-energy II(a;) = ij{k) with no extra 
MFL or SF frequency dependence. The parameters are presented in the text. 



14 



